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Abstract. Within the Nagel- Schreckenberg traffic flow modei we con- 
sider the transition from the free flow regime to the jammed regime. We 
introduce a method of analyzing the data which is based on the local 
density distribution. This analyzes allows us to determine the phase dia- 
gram and to examine the separation of the system into a coexisting free 
flow phase and a jammed phase above the transition. The investigation 
of the steady state structure factor yields that the decomposition in this 
phase coexistence regime is driven by density fluctuations, provided they 
exceed a critical wavelength. 

1 Introduction 

Over the past few years much attention has been devoted to the study of traffic 
flow. Since the seminal work of Lighthill and Whitham in the middle of the 50 's 

many attempts have been made to construct more and more sophisticated 
models which incorporate various phenomena occurring in real traffic (for an 
overview see ^|). Recently, a new class of models, based on the idea of cellular 
automata, has been proven to describe traffic dynamics in a very efficient way 
[0. Especially the transition from free flow to jammed traffic with increasing 
car density could be investigated very accurately. Nevertheless, besides various 
indications Q , no unique description for a dynamical transition has been found 
(see for instance [|]]6) and references therein). In this article we consider a method 
of analysis which allows us to identify the different phases of the system and to 
describe the phase transition, i.e., considering the fluctuations which drive the 
transition, and determining the phase diagram. 

We consider a one-dimensional cellular automaton of linear size L and N 
particles. Each particle is associated the integer values Vi € {0, 1, 2, w max } 
and di € {0, 1,2,3, ...}, representing the velocity and the distance to the next 
forward particle Hj. For each particle, the following update steps representing 
the acceleration, the slowing down, the noise, and the motion of the particles 
are done in parallel: (1) if Vi < di then Vi — > Min{i>i + l,u max }, (2) if Vi > di 
then Vi — ► di, (3) with probability P Vi — > Max{vi — 1,0}, and (4) r,i — > ^ + v i: 
where denotes the position of the z-th particle. 
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2 Simulation and Results 

Figure [l] shows a space-time plot of the system. Each dot corresponds to a 
particle at a given time step. The global density p g — N/L exceeds the critical 
density and jams occur. Traffic jams are characterized by a high local density of 
the particles and by a backward movement of shock waves . One can see from 
Fig. ^ that in the jammed regime the system is inhomogeneous, i.e., traffic jams 
with a high local density and free flow regions with a low local density coexist. 
In order to investigate this transition one has to take this inhomogeneity into 
account. 

Traditionally one determines the so-called fundamental diagram, i.e., the di- 
agram of the flow vs the density. The global flow is given by, <P = p g (v) , where 
(v) denotes the averaged velocity of the particles. This non-local measurements 
are not sensitive to the inhomogeneous character of the system, i.e., the infor- 
mation about the two different coexisting phases is lost. In the following we 
consider a method of analysis which is based on the measurement of the local 
density distribution p(p) || . The local density p is measured on a section of the 
system of size 8 according to 

1 N 

y i— 1 

The local density distribution p(p) is plotted for various values of the global 
density p g in Fig. g. In the case of small values of p g , see Fig. ||a, the particles 
can be considered as independent (see below) and the local density distribution 




Fig. 1. Space-time plot for w max = 5, P = 0.5, and p g > p c . Note the separation 
of the system in high and low density regions. 
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Fig. 3. The local density distribution p(p g ,p) as a function of the global den- 
sity (horizontal axis) and local density (vertical axis), respectively. The colors 
correspond to the values of the probability p{p g , p), increasing from black to red. 
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Fig. 4. The structure factor S(k) for P = 0.5, %ax = 5 and for various values 
of the density p. The dashed line marks the characteristic wavelength Ao of the 
free flow phase. 



is simply Gaussian with the mean values p g and a width which scales with \f&. 
Increasing the global density, jams occur and the distribution displays two dif- 
ferent peaks (Fig. ^c). The first peak corresponds to the density of free particles 
and in the phase coexistence regime the position of this peak does not depend on 
the global density (see the dashed lines in Fig. ^|). The second peak is located at 
larger densities and characterizes the jammed phase. With increasing density the 
second peak occurs in the vicinity of the critical density p c (Fig. ||b) and grows 
further (Fig. ||c) until it dominates the distribution in the sense that the first 
peak disappears (Fig. ||d). The two peak structure of the local density distribu- 
tion clearly reflects the coexistence of the free flow and jammed phase above the 
critical value p c . In Fig. || we present the probability distribution as function of 
the global and local density. Above a certain value of the global density p g the 
two peak structure occurs. The behavior of the first peak yields a criterion to 
determine the transition point || and one gets p c = 0.0695 ± 0.0007 for P = 0.5 
and w max , respectively. 

In order to describe the spatial decomposition of the coexisting phases we 
measured the steady state structure factor M 



S(k) 



ikr 



r=l 



(2) 



where r/(r) — 1 if the lattice site r is occupied and r](r) = otherwise. In Fig. || 
we plot the structure factor S(k) for the same values of the global density as in 
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Fig. 5. The phase diagram of the Nagel-Schreckenberg model. Note that in the 
non-deterministic region < P < 1 the density of the maximum flow exceeds 
the density of the transition point. 



Fig. |2], i.e., below, in the vicinity, above and far away of the transition point. It is 
remarkable that S(k) exhibits a maximum for all considered values of the global 
density at fco ~ 0.72 (dashed lines in Fig. ^). This value correspondence to the 
characteristic wave length Ao = of the density fluctuations in the free flow 
phase. The steady state structure factor is related to the Fourier transform of the 
real space density-density correlation function. The wave length Ao corresponds 
to a maximum of the correlation function, i.e., Ao describes the most likely 
distance of two particles in the free flow phase. For low densities the structure 
factor is almost independent of the density and displays a minimum for small 
k values indicating the lack of long-range correlations. Crossing the transition 
point the smallest mode S(k = ^) increases quickly. This suggests that the 
jammed phase is characterized by long-range correlations which decay in the 
limit p g » p c algebraically as one can see from the log-log plot in Fig. |^d. 

Up to now we only considered the case P — 0.5. The phase diagram in Fig. || 
shows the P dependence of the transition density p c . f denotes the free flow 
phase and f+j corresponds to the coexistence region where the system separates 
in the free flow and jammed phase. The dashed line displays the P dependence 
of the maximum flow obtained from an analysis of the fundamental diagram Q . 
The critical densities p c , where the phase transition takes place, are lower than 
the density values of the maximum How. Measurements of the relaxation time, 
which is expected to diverge at a transition point 0] , confirm this result |J (see 
Fig. |). But one has to mention that the determination of the critical density 
via relaxation times leads in the coexistence regime f+j to unphysical results, in 
the sense that the relaxation time becomes negative (|]|] ■ 
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3 Conclusions 

In conclusion we have studied numerically the Nagel-Schreckenberg traffic flow 
model using a local density analysis. Crossing the critical line of the system 
a phase transition takes place from a homogeneous regime (free flow phase) 
to an inhomogeneous regime which is characterized by a coexistence of two 
phases (free flow traffic and jammed traffic). The decomposition in the phase 
coexistence regime is driven by density fluctuations, provided they exceed a 
critical wavelength A c . 
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